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Abstract 

In this paper we review Penrose's Weyl curvature conjecture which 
states that the concept of gravitational entropy and the Weyl tensor is 
somehow linked, at least in a cosmological setting. We give a description 
of a certain entity constructed from the Weyl tensor, from the very early 
history of our universe until the present day. Inflation is an important 
mechanism in our early universe for homogenisation and isotropisation, 
and thus it must cause large effects upon the evolution of the gravitational 
entropy. Therefore the effects from inflationary fluids and a cosmological 
constant are studied in detail. 

1 The arrow of time and gravitational entropy 
in the context of a cosmology 

There is a strange omission in the traditional version of the second law of ther- 
modynamics (SLT). It does not take gravity into account. However the exis- 
tence of the arrow of time is usually explained with reference to SLT. And most 
discussions of the origin of the arrow of time appeal ultimately to the initial 
condition and the evolution of the universe. Gravity plays an essential role for 
this evolution. 

As pointed out by P.C.W. Davies |l], |2| there seems to be a paradox that 
the material contents of the universe began in a condition of thermodynamic 
equilibrium, whereas the universe today is far from equilibrium. Hence the 
thermodynamic entropy has been reduced in conflict with SLT which says that 
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the entropy of the universe is an increasing function of cosmic time. This is due 
to the tendency of self-gravitating systems irreversibly to grow inhomogeneous. 

In order to include this gravitational effect into a generalized version of SLT, 
one has to define a gravitational entropy. Several tentative definitions have been 
given based on Penrose's Weyl curvature hypothesis || 0, |j. 

The Weyl curvature tensor vanishes identically in the homogeneous and 
isotropic Friedmann-Robertson- Walker universe models. It was suggested by 
Penrose j| f|, [j| to use this tensor as a measure of inhomogeneities of the uni- 
verse models. The Weyl curvature scalar is expressed by scalars constructed 
from the Riemann curvature tensor as follows 

C a ^ s C aPlS = R^Ra^s - 2R aP Ra[3 + \R 2 (1) 

Penrose has given several formulations of the conjecture. A version that 
emphasizes the entropy aspect was given in ref. || p. 178. It seems that in 
some way the Weyl tensor gives a measure of the entropy in the space- time 
geometry. The initial curvature singularity would then be one with large Ricci 
tensor and vanishing Weyl tensor (zero entropy in the geometry); the final 
curvature singularity would have Weyl tensor much larger than Ricci tensor 
(large entropy in the geometry). 

In order to quantify Penrose's conjecture Wainwright and collaborators |], Q 
have suggested that the quantity 

R^R^u 1 ' 

may represent a "gravitational entropy" , at least in a cosmological context with 
a non-vanishing Ricci-tensor. 

We could of course ask ourself whether the entropic behavior of the Weyl 
curvature invariant (or an invariant composed thereof) results because the Weyl 
tensor is directly related to the gravitational entropy, or whether it is only a side- 
effect because the universe evolves towards a state of maximal entropy. Hence, 
even though the results may show that the Weyl curvature invariant has an 
entropic behaviour, it is by no means a proof that the Weyl tensor should be 
identified as the "gravitational entropy" . 



2 The behaviour of the Weyl tensor 

In an earlier article Q we discussed the Weyl curvature tensor in a homogeneous 
but anisotropic model as well as in an inhomogeneous model. Our survey will 
be summarized in this and the next sections, as well as bringing new arguments 
into the discussion. 

As a basis for our study, we used two different models which induce Weyl 
curvature effects in two conceptually different ways. One was the anisotropic 
Bianchi type I model, the other was the inhomogeneous Lemaitre-Tolman model. 
The Weyl curvature conjecture has been investigated in the Szekeres cosmolog- 
ical model, that generalize the Lemaitre-Tolman model, by W.B. Bonnor ||. 
Inhomogeneous modes are in general local modes for general relativity while 
anisotropic modes are global modes. The global topology and geometry of our 
universe has significant consequences for the possibilities for anisotropic modes 
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00, (f§ The role of inhomo aeneous modes can be considered as more 
local in origin, but the effect on the global geometry and topology are not known 
in detail. Hence, it is important to study both an inhomogeneous model as well 
as a homogeneous model to get a more complete picture of the behaviour of the 
Weyl tensor in generic cosmological models. 



2.1 The Bianchi type I model 

The Bianchi type I model is the simplest anisotropic generalisation of the FRW 
models. It has flat spatial sections, and hence, can be a good candidate for the 
universe we live in. The metric for this model can be written as 

ds 2 = ~dt 2 + e 2a [e 213 ] dx'dx 1 (3) 

where /3 = diag(/3+ + %/3/3_, (3+ - V3/3-, -2/3+). In Q we solved the Einstein 
field equations for dust and a cosmological constant A. By introducing the 
volume element v = e 3a , the Einstein field equations gives the following equation 
for v: 

v 2 = 3Av 2 + 3Mv + A 2 (4) 

where M is the total mass of the dust, and A is an anisotropy parameter. The 
equations for fj± are in terms of the volume element 

4 = f (5) 

where the constants a± are related to A via a\ +a 2 _ = A 2 . It is useful to define 
an angular variable 7 by a+ = Asin(7 — it/ 6) and a_ = Acos(7 — tt/6). An 
important special case of the Bianchi type I is the Kasner vacuum solutions. 
The Kasner solutions are characterised only by the angular variable 7. These 
solutions have a Weyl scalar 

Ifi A 4 

(CPTCaw), = --f (1 - 2zcos3 7 + z 2 ) (6) 

where z = v/A and z = 1 for the Kasner solutions which have vanishing Ricci 
tensor. So in a sense, the Kasner solutions are the anisotropic counterpart to 
the inhomogeneous Schwarzschild solution. 

The Weyl tensor will decay as the volume expands, even at late times. One 
should expect that this entity would increase if it represents gravitational en- 
tropy. But is decreases monotonically, and hence, it is doubtful that it is the 
correct measure in these models. 

Inserting dust and a cosmological constant will not have a significant effect 
on this decreasing behaviour. If M is the mass of the dust inside the volume v, 
then eq. (H) still holds but with 



3Mv 3Av 2 
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Hence, except for the special case]] 7 = 0, the Weyl tensor diverges as v 4 as 
v — ► 0. 

For the Bianchi type I universe models the entity P 2 defined in eq. (0), 
turns out to be 

2 _ A A A l + z 2 -2zcos3 7 ( ] 

{ H 27 v 2 M 2 + 2AMv + 4A 2 v 2 ' [ ' 

Also this measure for the Weyl entropy diverges as v — > 0. 
2.2 The Lemaitre-Tolman models 

Let us now consider the inhomogeneous Lemaitre-Tolman (LT) models. The 
line element for the LT models can be written as 

ds 2 = -dt 2 + Q 2 dr 2 + R 2 (dd 2 + sin 2 9d4> 2 ) (9) 

where Q = Q(r, t) and R = i?(r, t). In jl5| we used this model and investigated 
the solutions of the Einstein field equations where the spacetime contains a 
cosmological constant A and dust. The equations then turn into 

R' = FQ 

where F = F(r) is an arbitrary function, and 

^RR 2 + ^(l-F 2 )R-^R 3 =m. (10) 

The function m = m(r) is given by the integral 

m{r) = [ A-K P R 2 R'dr, (11) 
Jo 

where prime denotes derivative with respect to r and p is the dust density. 
Hence, m(r) can be interpreted as the total mass of the dust inside the spherical 
shell of coordinate radius r. Inverting the equation for m(r), the dust density 
can be written in terms of R and m(r): 

m' 

Ann = (12) 
H R?R' y ' 

It is also useful to define the mean dust density function p(r) by the relation 

m(r) = ^npR 3 . (13) 

Interestingly, the Weyl curvature scalar can now be written quite elegantly as 

1fi 2 

(C a ^ s C afflS ) LT = -fn 2 (p- P ) 2 . (14) 



lr The special case of 7 = is somewhat interesting. For 7 = it can be shown that the 
Kasner solution is just a special part of Minkowski spacetime, which has no singularities. It 
seems a bit odd that for any 7 / we will have an initial singularity, while for 7 = we 
do not have one. However, if we compactify the spatial sections this oddity disappears. The 
model has a singularity at v = even for the case 7 = 0! 
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Ever expanding solutions 



Figure 1: Level curves of E(R, R) = \R 2 - ^ - \R 2 



This relation provides us with a physical interpretation of the Weyl scalar in 
the LT models. It is just the difference between the mean dust density and the 
actual dust density. Also, the Weyl tensor is everywhere zero, if and only if 
p = p. In this case the LT models turn into the FRW models with homogeneous 
spatial sections. The FRW models are conformally flat, hence they have a zero 
Weyl tensor. 

Let us return to the equations of motion for the LT model. The mass equa- 
tion (|l0| ) can be written as an "energy" equation: 

l -R 2 + V{ ri R) = E{r) (15) 
where the "potential" V and the "energy" E are given by 

R 6 

E = -\{l-F 2 ) (16) 

This energy equation may be integrated and solved exactly. A summary of the 
results and a qualitative description of the physical meaning of the solutions 
is given in pq| . The solutions may be written in terms of the Weierstrass' 
elliptic functions |l7| . The actual expressions are not very informative unless 
the reader has massive knowledge of these elliptic functions. However a lot of 
qualitative information can be extracted from simple classical considerations. 

The classical solutions will move on level curves of the energy function 
E(R,R) = ii? 2 — ^ — ^R 2 , since the total energy E is independent of t. 
In figure [l] the level curves of a typical energy function are drawn. If A > 
there will exist a saddle-point of the energy function. This saddle point will 

be at R — (^p) 3 ,R — where the energy function will have the value 
E s = — i(9m 2 A)^. The saddle point solution is static, and is as a matter 
of fact the Einstein static universe^. If E < E s , the solutions fall into two 
distinct classes: 



2 This is easily seen if we define the new radial variable to be _R. This can be done since 
we have to assume that m'(r) > on physical grounds. In the case m'(r) = 0, the metric 
become degenerate. 
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1. Schwarzschild-like solutions: These solutions expands, but they do not 
possess enough energy to escape the gravitational collapse, so they end as 
black holes. If m is constant these solutions are those of a Schwarzschild 
black hole in Lemaitre coordinates. 

2. de Sitter-like solutions: solutions where the universe evolves approxi- 
mately as that of de Sitter solutions with positively curved hypersurfaces. 

If E > E s the (test) matter has enough energy to escape the gravitational 
collapse (expanding solutions) or enough energy to prevent the gravitational 
repulsion from the cosmological constant (contracting solutions). 

Near the initial singularity, both p and p diverge. Unless they are identically 
equal, the Weyl tensor will diverge near the initial singularity. Explicitly we 
have near the initial singularity 

^(p- p ) = ^ 2 -$ (17) 

where we have used that near the initial singularity, we can approximate the 
solutions with 

fl» Qm) 3 (t-io(r))t. (18) 

The free function io( r ) is the big bang time. To avoid intersecting world- lines 
we have to assume t' < 0|p8| , thus (p — p) > 0. The Weyl scalar will diverge as 
R — ► unless t' — 0. Also the entity P will diverge in general near the initial 
singularity. P is given by 



\ 47T/3 J \ A-Kp 

which near the initial singularity can be approximated by 

Hence, P diverges as t — ► t unless t' — 0, i.e. unless the big bang is homoge- 
neous. 

This does not prove to be a very promising behaviour for the WCC. The 
curvature scalars diverge near the initial singularity, and hence, come in conflict 
with the WCC. But let us analyse the situation in the LT a bit more carefully. 
As the universe expands, both p and p will decrease. They both decrease from 
an infinite value at the initial singularity. As the universe expands the value of 
p decrease as t~ 2 while p decreases as t" 1 close to the initial singularity. Hence, 
the value of p approaches the value of p by a factor or t , but since they both 
diverge, the Weyl tensor and P diverge as well. 

However, even though R. Penrose and S. Hawking [[l9) showed that according 
to the general theory of relativity the big bang must have started in a singularity, 
it has later been emphasized that the physical universe must obey not only the 
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Figure 2: The plot of the function \P\ = T~ 
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general relativistic laws of nature, but also the quantum mechanical laws. Hence, 
the initial singularity is a fiction which not corresponds to physical reality. The 
classical laws can be applied only after the Planck time. We should therefore 
study the behaviour of P 2 not only in the limit t — > 0, but rather at a very 
small cosmic time. Choosing the origin of cosmic time T at t = to we introduce 
T = t - i (r). Hence close to T = 



\P\ = 2- 



J*o[ 
T 



(21) 



which implies that [d\P\/dT] T >0 < 0. Bonnor considered the model with 

1 — F 2 > and a vanishing cosmological constant, and found the opposite 
result. However he restricted his investigation to models with homogeneous 
initial singularity, i.e. t' = 0. As we will show this is an exceptional case. 
The behaviour of more general models with t' ^ is different. To see what is 
actually happening if we push t' towards zero, we utilize that close to T = 
the models with 1 — F 2 > and A = behave according to 



R 3 = ^mT 2 (l-e(r)T^ 



(22) 



where e(r) is assumed to be a small function of r. From eqs. (|l2"|), ( p"3| ) and ( p^ ) 
with A = we get 



IPI 



2 ( mB! 1 



V3 V m ' R 



(23) 



Using (E2h we obtain 



m 



2m\t'\ T 



leT- 
3 



e'Tf 



(24) 



This function is plotted in fig. 

Note that this function diverges for T 



0, obtains a minimum for some 



small T and increases thereafter, 
goes towards T = 0. 



As in is driven towards zero this minimum 
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3 The Gravitational entropy revisited 



Let us consider a finite 3- volume V comoving with the cosmic gas in our space 
time. According to the first law of thermodynamics the matter entropy Sm and 
the internal energy U will evolve as: 

TdS M = dU+pdV. (25) 

If the matter content is dust then p = 0. There are a couple of things to 
note. Firstly in a dense dust cloud, we expect the internal energy of the dust 
to be large, thus we expect the entropy to be large. Secondly, the entropy is 
an increasing function of the volume. Let us therefore consider a co-moving 
volume V in our spacetime. What could the expression for a gravitational 
entropy be? From the previous sections we noticed that an LT model with 
p = p is homogeneous. The quantity is an inhomogeneity measure in the 
LT models. In the absence of a cosmological constant we notice that 

The sign of P is here chosen so that the configuration p < p is associated with 
P < while the more realistic configuration p > p has P > corresponding to 
eq. (^6|) positive. Let us consider the entity defined by: 

S= I PdV. (27) 
Jv 

Introducing co-moving coordinates x l we write dV = ^/hd 3 x where \fh is the 
3- volume element. The integration range is now constant as a function of time 
and if we integrate over a unit coordinate volume which is so small that the 
integrand is approximately constant, we may write 

S= I PdVfzPVh (28) 
Jv 

This formulation of the gravitational entropy is an intuitive one, but it is not 
a covariant formulation. To find a covariant formulation we define the entropy 
current vector by Q 

* = ,su + ip (29) 

where u and y are orthogonal, s is the entropy density, u is the material flow 
vector (uu = —1) and ip is the entropy flux. The second law of thermodynamics 
can now be expressed as 

> 0. (30) 
Writing this in a local coordinate system, the divergence is 



l=d, (vIsH • (31) 



In our case we take the entropy density proportional to P and assume a vanishing 
entropy flux. Hence, 

* cx Pu. (32) 
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In comoving coordinates u = e t so that 

ex l=d t (VhP) = -j=S. (33) 

For the "second law to hold" , we have to check whether the entity S = \[hP 
is increasing. If we take the corresponding one-form of the entropy vector, and 
take the dual of this form (Hodge dual), we obtain the entropy three-form which 
is given by the contraction of the space-time volume form e with the entropy 
current current vector 

ui = type. (34) 

In our case, this three-form is simply 

oj = VhPdx A dy A dz (35) 

Hence, if the coordinates are comoving, then it is the component of this three- 
form that is increasing as long at the matter and fields obeys the SEC. This 
explains why the entropy should scale as the volume. It is usually this entropy 
three-form we are thinking of. It has a more intuitive behaviour than the entropy 
current vector. 

We will therefore use S in the further to study the WCC. The square of 
the Weyl tensor and the entity P do not, as we have seen, capture the entropic 
behaviour properly. They both diverge at the initial singularity and thus, cannot 
be a proper measure of the gravitational entropy. The entity iS on the other 
hand, seems to be a more promising candidate for the Weyl entropy. 



4 The behaviour of S 

4.1 The Lemaitre-Tolman model 

Let us start with the LT model. The motivation for studying the entity S given 
in came from the LT models. Therefore it is natural to start with this 
model. In all the LT models, we can approximate the behaviour near the initial 
singularity, R — ► with 

S LT = 2V3^(m'(t-t )-2t' ). (36) 

For physically realistic spacetimes, t' < 0, so as t — > t' , S is positive and finite. 
Note also that S is increasing, as suggested by the WCC. In the absence of a 
cosmological constant and if F = 1, this will be the exact expression for S. If a 
cosmological constant is present, the universe will eventually go into a de Sitter 
phase if the universe is allowed to expand for ever. 

For the sake of illustration, it is useful to investigate a specific case. Let us 
choose F 2 = 1. The energy equation eq. (tt5h, can now be solved to yield 



\H{t-t Q {r)) 



(37) 



where H = yj A/3. In general, the late time behaviour of the LT model (if the 
universe has one) is that of a de Sitter universe. This exact solution provides 
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Figure 3: The evolution of S in the LT models is dependent on the value of A. 
The upper graph has a small value of A which makes S increase to a large value. 
The lowest graph has a large value of A. 



us with a solution which connects the initial singularity with the late time de 
Sitter era. The initial behaviour of R is the same as in eq. (|l8|). The late time 
behaviour of R is 
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R=[^Ye H ^\ (38) 



2A 

Hence, at late times the entity S is approximately constant: 

S= 2 ^. (39) 

Aa 

At late times, S approaches a constant value which is inversely proportional to 
H. The larger the value of the cosmological constant is, the smaller the value 
of the final value of S. In fig. |^ we have plotted, using the exact solution (p7|), 
the evolution of S for three different values of A. 

The volume under consideration expands at late times exponentially as V oc 
Vo(r) exp (3Ht), where H is the Hubble parameter. Hence, the entropy per unit 
volume will decrease exponentially: 

£ oc e-™\ (40) 

The entropy in a unit volume will therefore decrease rapidly during a de Sitter 
stage. This is a very important consequence of inflation. At the exit of the 
inflationary period, the entropy in a unit volume is very low, there are very 
little inhomogeneities left of the primordial ones. 

Let us summarize the evolution of the entity iS in the LT models: 
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• Large A and ever expanding: In the initial epoch the dust dominates 
and S is increasing linearly in t. The universe is becoming more and more 
inhomogeneous. After the universe has grown considerably, the cosmolog- 
ical constant becomes dominant, S stops growing and if the cosmological 
constant is large enough, it evolves asymptotically towards a constant 
value. The universe is smoothened out. This scenario corresponds ap- 
proximately to the lowest two graphs in Fig. ^. 

• Small A and ever expanding: Again the dust dominates initially. The 
cosmological constant is too small to make S decreasing. The entity S is 
ever increasing but is bounded from above by a relatively large constant 
value. 

• Zero A and ever expanding: The S will again be ever increasing and 
will asymptotically approach a function f(t) = c + bt p where c and b are 
constants and p = 3 iff F 2 > 1 and p = 1 iff F 2 = 1. 

• Recollapsing universe: Due to the dust term, the final singularity will 
not be similar to the initial singularity. Hence this entity is asymmetric 
in time for a recollapsing universe. 

In the LT models we see that S behaves in agreement with the WCC. 

Let us investigate more carefully the Schwarzschild spacetime, where the 
whole motivation of gravitational entropy comes from |2^, |2^] . The Schwarzschild 
spacetime is a special case of the general LT models; it has m(r) = constant 
and a vanishing Ricci tensor. If we look at the entity iS in the region outside 
the Schwarzschild singularity, S will diverge. This is in some sense the maxi- 
mal possible value of S, the Weyl tensor is as large as possible and the Ricci 
tensor is the smallest as possible. Thus at this classical level it seems that 
the Schwarzschild spacetime has the largest possible S which is a good thing 
if one wants to connect iS with the entropy of the gravitational field. Outside 
a real black hole, R^ u — is probably impossible. Even though the classical 
vacuum has vanishing Ricci tensor, the quantum vacuum will probably not have 
R^v = 0. The quantum fields will fluctuate and cause the expectation value of 
the square of the Ricci tensor to be non-zero: (R^R^) ^ 0. Hence, there will 
probably be an upper bound of how large S can be, even in a vacuum. 

4.2 The Bianchi type I model 

In the Bianchi type I with dust, the explicit expression for the entity S is 

rr 2 A 2 ( l + z 2 -2zcos37 V 
Sl = VhP= ^{M 2 + 2AMv + 4A 2 v 2 ) (41) 

where \[h = v. In a neighbourhood of the singularity we can make a Taylor 
expansion to first order in v: 



„ 9 8 A 4 , J ( 3M 2A 



We interpret m = M/A as the dust density in coordinate space || . Hence, 
as we go towards the initial singularity S will be finite. For the vacuum case, 
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m = and A = 0, and the solutions are the Kasner solutions. The Kasner 
solutions have a vanishing Ricci tensor, but a non-zero Weyl tensor. Hence, as 
in the Schwarzschild case, the entity S will diverge. 

From the Taylor expansion we see that, if 3to 2 > 4 A then S will increase 
immediately after the initial singularity. If 3to 2 < 4A then S will decrease. 
This can be understood as follows. If the cosmological constant is to large, 
the universe will increase too rapidly initially to allow the dust to contribute 
significantly to the anisotropy. Initially the universe will become more and more 
isotropic. In the presence of a cosmological constant, the universe will eventually 
enter a de Sitter phase. Again S will evolve towards a constant value given by 

lim Si = -. 

v ^°° 3A2 

If the cosmological constant vanishes S is monotonically increasing: 

2A ( 9 \ ^ 

Si = [{2 + 3rot)(l - cos3 7 ) + -m 2 t 2 J 

Note that for large t we have Si = -j=t, which is according to the same power 

law as the F 2 = 1 and A = case of the LT model (compare with eq. (|3q)). 

To summarize our investigation of S in the Bianchi type I model we can 
say that the entity S at the initial singularity is a constant determined by the 
inverse of the dust density. For 3m 2 > 4A it will increase immediately after 
the initial singularity. In most cosmological considerations it is assumed that 
the cosmological constant is small. The exception is in the inflationary era in 
which the vacuum energy dominates over all other matter degrees of freedom. 
The inflationary era will smooth out anisotropics as well as inhomogeneities, 
and the behavior of S in this case is therefore expected. In the de Sitter limit 
S will asymptotically move towards a constant. Large A means small value, 
while small A corresponds to a large value. This is in full agreement with the 
LT models. It is also interesting that the entity S is very sensitive to different 
matter configurations. This makes it a lot easier to check whether the S has 
the right entropic behaviour. 

The question now arises: How generic is this behaviour? Does the entity 
behave in the correct way for all physically realistic models? 

We know that our universe today is close to homogeneous on a scale larger 
than a billion light years. As mentioned earlier, the FRW models are conformally 
flat and hence, they have a zero Weyl tensor. So a question would be, at late 
times when the universe is close to isotropic and homogeneous, does the entity 
S still behave in the correct manner? Does it still increase? We saw that at late 
times, both the LT model and the Bianchi type I model, did behave correctly 
for A = even though they both isotropise and evolve towards homogeneity. 

Let us choose a more general model, a flat model that allows for inhomo- 
geneities and anisotropies. We will assume that the model at late times asymp- 
totically evolve towards a FRW model. In |24|, ^5), Barrow and Maartens in- 
vestigated the general equations of motion for such models with A = 0. In the 
velocity dependent regime they derived some approximate solutions to the field 
equation close to a FRW model, having both inhomogeneities and anisotropy. 
Also an anisotropic stress of the form = \^p r where A Ml/ is a constant 
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matrix, was included. Their result was that the late time behaviour of such a 
model is consistent with a perfect fluid where the equation of state parameter 7 
for the perfect fluid obeys 7 < 4/3. For 7 > 4/3 the FRW model is unstable at 
late times. In the latter case one can get for instance different stable anisotropic 
magnetic solutions |^6|, ^?J. 

Let us therefore assume that 2/3 < 7 < 4/3 (so that the matter obeys the 
SEC). The ratio of the Weyl tensor squared and the Ricci tensor squared is in 
this case 

p2 = ^ aM& ^ C ( t -2n) + 0( t -3») + {t- An ) (42) 

ft ftfiu 

where n — 2 ^ 4 3 ^ 37 - > for n =/= 4/3. In the radiation case, 7 = 4/3, we get a 
logarithmic decay of P 2 : 

P ^/3«(InV ^ 

Thus in the range of validity of the assumption 7 < | , this entity will decrease 
in the future. However, the entity S shows more promising behaviour 

S = PVh = a 3 P oct^" n (44) 

for 2/3 < 7 < 4/3, while the 7 = 4/3 case yields 

S oc (lni)- 1 ^. (45) 

Hence, S increases as long as 

§<7<|- (46) 

This entity increases in the future (as any entity describing entropy should do) 
and tells us that even if the universe itself asymptotically goes towards isotropy, 
the entropy of the gravitational field actually increases (if we should believe the 
WCC). 

Note that this result is slightly different than the case where no anisotropic 
stress is present. If the anisotropic stress is not present, we get S oc t at late 
times for all 7. 

We should also mention a work done by Hervik [p8| which investigates the 
evolution of the Weyl curvature invariant for generic solutions of spatially homo- 
geneous models containing a 7-law perfect fluid with 7 > 2/3. The conclusion 
was that all spatially homogeneous models, except for sets of measure zero, had 
an increasing S at late times. 



4.3 Inflation 

As long as the matter obeys the strong energy condition (SEC), matter will 
behave more or less attractive. Hence, when the SEC is fulfilled, we should 
expect the gravitational entropy to increase. 

During inflation the SEC is violated and gravity is not necessarily attractive. 
We have already seen how the inclusion of a cosmological constant could alter the 
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behaviour of S. The late time behaviour of S in the presence of a cosmological 
constant is that S evolves approximately as a constant. The constant itself is a 
decreasing function of A (see eq. (|39|)). 

What happens for a more general inflationary fluid? Consider a perfect fluid 
with equation of state 

V = (7 ~ 1)P (47) 

for < 7 < |. This perfect fluid will violate the strong energy condition and 
will in general cause a power law inflation. Note that the case 7 = can be 
considered the same as including a cosmological constant. In the FRW cases, 
all the models will now have a late time behaviour similar to the flat case. The 
FRW models will become inflationary and will be dominated by this fluid at 
some stage in the future (provided that the universe is ever-expanding). To 
simplify, we will therefore consider a flat universe, and perturb the FRW flat 
universe model with 7 < 2/3. Following |24|, [25), the anisotropic stresses will not 
dominate the shear modes at late times. The shear will under these assumptions 
decrease as 

<7 M „ OC t~v (48) 

while the Hubble parameter is the same as in the FRW case (to lowest order) 

H = i- t -i. (49) 
37 

The late time behaviour of P 2 is now 

p2 =Q ( 7 2 t -*^) (50) 

which is decreasing for all < 7 < 2/3. The late time behaviour of S is 

S = VhP = O (-ft) (51) 

which is increasing linearly in t at late times. For 7 = we recover the cosmo- 
logical constant case where S is approximately a constant at late times. Note 
also that this result coincide with cq. ( ]44| ) for 7 = 2/3. 

Hence, during the inflationary period, S increases much slower than for 
ordinary matter ( S cx 1 3 for dust and S oc (lni) -1 ^ for radiation). 



5 Using Quantum Cosmology to determine the 
initial state 

Maybe the closest to a quantum theory of gravitation that have been obtained, 
is what we call quantum cosmology (QC). We will in this section show how one 
might be able to determine the likelihood of a certain initial state to occur. QC 
is perhaps best described as a "theory of initial conditions" . We will again use 
the results of the previous paper to try determine which of the initial states are 
more probable. 
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In QC, the wave function of the universe satisfies the Wheeler-DeWitt (WD) 
equation: 

(- G ^^ +y ^* = ° < 52 > 

where Gijki is called DeWitt's supermetric, and hij is the metric on the 3- 
dimensional spatial hypersurfaces. The potential term V(hij) consists of the 
Ricci scalar of the three-dimensional hypersurfaces and possibly matter poten- 
tials and a cosmological constant. 

We want now to calculate the expectation value for the Weyl scalar, in 
these models. All of the curvature invariants have to go over to their respective 
curvature operators. Especially, the entity S goes over to the curvature operator 
S. The question is now, what is more likely, creation of a universe with large 
expectation value of S or a universe with a low value? 



5.1 The LT model 

In the paper |15j we considered semi-classical tunneling wave functions which 
are solutions of the WD equation. The universe was tunneling from a matter- 
dominated universe classically confined to a finite size, into a de Sitter like 
universe. After tunneling across the classically forbidden region the universe 
became A-dominated, similarly to an inflationary model of the universe. Thus 
we have a "flow" from matter-dominated universes towards de Sitter like uni- 
verses. 

The actual expectation values of entities like C a ^ lS C a p 1 s 1 P 2 and S are not 
explicitly obtained for these models, because the actual calculations suffer from 
"endless" expressions and highly time-consuming quantities. We will therefore 
give a more general description of the evolution of the Weyl tensor for the LT 
models. 

It would be useful first to recapitulate some of the discussion done in |il], || . 
First of all we discussed tunneling wave functions in the WKB approximation. 
In the WKB approximation we assume that the wave function has the form 
^wkb — exp(±iS), where S will to the lowest order satisfy the Hamilton- 
Jacobi equation: 



SS \ 2 F' 2 



SR F 4 



2mR-R 2 {l-F 2 ) + ^R 4 



= (53) 



Here S/8R denotes the functional derivative with respect to the function R. 
If we assume that S = J o~(r)dr, the resulting equation will be the Hamilton- 
Jacobi equation for a point particle with action <r(r) (r is only a parameter 
and the functional derivatives turn into ordinary partial derivatives). In the 
Hamilton- Jacobi equation the functional S turns out to be the action at the 
classical level. Since the classical action can be written as an integral over r 
the assumption S = j a{r)dr is therefore reasonable at the lowest order WKB 
level. We can interpret the action a as the action of a point particle moving in 
a potential V(R) = ^ [-mR+ |(1 - F 2 )R 2 - jrR 4 ] with zero energy. The 
WKB wave function ipwKB for the point particle can then be written ipwKB — 
exp(±zcr). The two WKB wave functions can therefore be related by ^wkb — 
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cxp(J dr In ipwKB ) ■ Finding first the wave function ip we can then relate its 
WKB approximation to ^wkb through ^wkb = exp(J dr In ipWKB ) • 

Let us now ask the question: Is it more likely for a universe with small Weyl 
tensor to tunnel through the classical barrier than a universe with a large Weyl 
tensor? The question is difficult to answer in general but we shall make some 
simple considerations in order to shed some light upon it. 

We assume that the dust density near the origin of the coordinates is larger 
than further out. We define the homogeneous mass function for a closed uni- 
verse (k = 1) as mh(r) = ^nphr 3 where ph is a constant. The constant ph is 
determined by demanding mh{r m ax) = ni(r max ). If the dust density is larger 
near the origin of the r-coordinate than for larger values of r then m(r) > mh(r). 
This will not in general change the size of the universe so we can look at the 
effects from m(r) alone. Since m(r) is greater in general for an inhomogcneous 
universe than for a homogeneous universe, we see that the potential barrier will 
be smaller for an inhomogcneous universe than for a homogeneous universe. 
Thus an inhomogcneous universe will tunnel more easily through the classical 
barrier than the homogeneous universe. Since an inhomogeneous universe will 
have a larger Weyl tensor than an almost homogeneous one, we see that uni- 
verses with large Weyl tensor tunnel more easily than those with a small Weyl 
tensor. 

If we look at the tunneling amplitude concerning effects from the A term, it 
is evident that larger A will yield a larger tunneling probability. In the initial 
era inhomogencitics will increase the value of S. We saw that an inhomogcneous 
state will tunnel more easily through the potential barrier than a homogeneous 
state. The largest tunneling probability amplitude thus occurs for universes with 
a large cosmological constant and large local inhomogeneities. From a classical 
point of view the value of S initially was large (but increasing thereafter), but 
as the universe entered the inflationary era the cosmological constant had a 
large value, hence the value of S at the end of the inflationary era was relatively 
small. 

In the initial epoch the universe is not believed to be dust dominated. The 
dust does not exert any pressure and dust particles do therefore not interact 
with each other. A more probable matter content is matter which has internal 
pressure. Even though gravitation tends to make the space inhomogeneous, 
internal pressure from the matter will try to homogenise the space. Since dust 
is the only matter source in our model, the model only indicates the tendency 
for gravitation itself to create inhomogencitics. 

Since the universes tunnel into a de Sitter-like state, the cosmological con- 
stant will rapidly dominate the evolution. The larger the cosmological constant 
the lower will the entity S be after the inflationary era ends. 



5.2 The Bianchi type I model 

We can write the general solution of the WD-equation for the Bianchi type I 
models as: 



*(«,/?) = / d 2 k [C(k)^(v)p(k) ( 



i/3-k 



(54) 



where ipz(v) = v 2 Wl,h(2Hv) is a particular solution of the WD equation (with 
dust), C(k) is a "normalizing constant", and p(k) is a distribution function in 
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momentum space. This distribution function satisfies the equation: 

J d 2 k\p(k)\ 2 = 1 

The function Wl^(x) is the Whittaker function. 

There is still a factoring problem in turning the classical entities to operators. 
Let us first investigate the expectation value of C al} ^ 5 C a p~ l s- In the expression 
for the Weyl tensor squared, there is a term 

^-Apl+P-f (55) 

where p± is the conjugated momenta to f3±. Upon quantization, the above term 
is replaced with the operator: 




In the expectation value of the square of the Weyl tensor, the above term will 
contribute with 

/^/ra*£(^+<fe) 2 tt _ i Jd 2 k\C(k)\ 2 \^(v)\ 2 \p(k)\ 2 \k\ 4 
fd 2 f3^*<$ ~ « 4 J d 2 k\C{k)\ 2 \^(v)\ 2 \p(k)\ 2 

Unless we have a delta- function distribution at k = 0; p(k) = S 2 (k), the contri- 
bution from this term to the Weyl curvature invariant will diverge as v~ 4 for 
small v. Thus, we have to conclude that, in the small v limit the expectation 
value of C al3lS C a f3 7 s goes as: 

{C a ^ 5 C aPlS ) cx 1 (57) 

just as in the classical case. 

Investigating the invariant R^R^, wc notice that things arc not so easy. 
The Ricci square also has a term which presumably would contribute with a 
% term. However, looking at the classical expression we see that the Ricci 
square is independent of the anisotropy parameter A. This indicates that at 
the classical level all terms involving the anisotropy parameters, have to cancel 
exactly. This is not the case quantum mechanically. In the quantum case 
operators do not necessarily commute. Hence there may be contributions from 
terms which classically would cancel each other. In other words, the fact that 
the classical vacuum has R^ v = 0, does not mean that the quantum vacuum 
has R^u = 0. 

We assume that (£j ) is a set of factor-ordering parameters which represents 
the "true" quantum mechanical system in such a way that £j = represents the 
classical system. With this parameterization of the factor ordering we would 
expect the Ricci square expectation value to be: 

(BTR^) = 4A 2 + 2A^ + ^ + ^ + 0(g) 
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where fi is some function of v which has the property: v s» 0, fj(v) ~ 
constant. Thus for small v and £j the expectation value would behave as 



and 



■ constant (58) 



The Weyl square divided by the Ricci square is in general finite as v — > for a 
quantum system. In some sense, the quantum mechanical effects renormalizcs 
the infinity that the classical system possesses at v = 0. The expectation value 
at v = is, however, strongly dependent on the factor-ordering. As the factor 
ordering parameters approach zero, the value will diverge. Quantum effects in 
the early epoch are essential for the behaviour of this entity near the initial 
singularity. As v — ► we expect the quantum effect to be considerable, thus 
expecting the factor-ordering parameters £j to be large. 

As indicated in the above discussion, the quantum mechanical expectation 
value of P will be lower and presumably finite at the initial singularity. Therefore 
the expectation value of S is also presumed to be considerably lower in the initial 
stages than its classical counterpart. 

Comparing different tunnelling amplitudes in the Bianchi type I model is dif- 
ficult and more speculative because the Bianchi type I universe has no classically 
forbidden region for A > 0. This causes the lowest order WKB approximation 
to be purely oscillatory. The lowest order WKB wave function will therefore 
be approximately constant. In the paper JL4| we did however construct under 
some assumptions a wave function which clearly peaked at small values of the 
anisotropy parameter. Thus these wave functions predicts universes that have 
a relatively low value of S. 



6 Cosmic evolution of the Weyl Curvature 

Let us now recapitulate how the evolution in the context of the WCC might 
have been. 

The universe was created in a rather arbitrary state. As the time ticked 
past the Planck time 10 -43 s a rather inhomogeneous universe appeared. QC 
suggests that this state was rather inhomogeneous, but it is more or less a guess 
how inhomogeneous the universe was at that state. Nevertheless, as the universe 
grew larger, the expectation value of the Weyl entropy increased initially. At 
what rate the Weyl entropy increased is highly uncertain, it depends very much 
on the true nature of our universe. It depends on what matter fields that were 
present, the topology of the universe, whether it was anisotropic or not, quantum 
effects etc. 

Nevertheless, at some time very short after the Big Bang, an enormous effec- 
tive cosmological constant appeared. The universe was driven unconditionally 
into an inflationary period. The Weyl entropy stopped increasing and began 
instead to evolve asymptotically towards a constant. Whether or not the Weyl 
entropy had earlier a higher value than at the exit of the inflationary regime, 
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time 



Figure 4: A schematic graph of the evolution of the Weyl entropy. 

is difficult to say. During the inflationary period, the universe was more or less 
in an adiabatic expanding state, the Weyl entropy was more or less constant. 
Since the universe was expanding exponentially during this period, the entropy 
per unit volume dropped exponentially. If the scale factor increased by a factor 
of 60 e-foldings during the inflationary epoch, then the entropy per unit volume 
of space would have decreased by a factor of 

- ~ e- 3 ' 60 « 10 - 78 . (59) 
so 

This is quite a drastic decrease, and shows how powerful inflation is when it 
comes to smoothing out the inhomogeneities and anisotropies of our universe. 

At the exit of the inflationary regime, the Weyl entropy had dropped by 
an enormous factor compared to what is would have been if no inflation had 
occurred. The Weyl entropy was very small, compared to the maximal possible 
value. The universe was more or less uniform and homogeneous. This homo- 
geneity can be seen in the CMB radiation today. However, small fluctuations 
in the spectrum can also be seen, reflecting the state of the universe 300 000 
years after the Big Bang. After the inflation, small seeds of inhomogeneities 
from the quantum fluctuations of quantum fields were the only thing left of the 
initial inhomogeneities. Nevertheless, these seeds were large enough to gradu- 
ally clump together and form galaxies and stars. The Weyl entropy began to 
increase again after the inflationary regime was over. 

The radiation-, and later the matter-dominated universe caused the Weyl 
entropy to grow steadily and firmly for almost 15 billion years. During the 
radiation era, there were only small inhomogeneities and anisotropies left of 
the primordial fluctuations. Using the calculations from the earlier section, the 
Weyl entropy increases as S oc (In t)~ 1 ti in the radiation era. 

At about t =10 000 years the radiation become sub-dominant. The universe 
evolved into a matter dominated era and the Weyl entropy increased as S oc 

4 

frs. Today it is still growing. Recent observations suggest that the universe 



19 



has entered a new era with accelerated expansion. The surprising fact that 
the universe appears to be in an accelerating state today, can be explained 
with the presence of a vacuum energy. If this is true, the Weyl entropy will 
increase steadily and asymptotically towards a constant. However, the late 
time behaviour of our universe is still quite speculative, and whether or not this 
vacuum dominated period will persist, is very uncertain. It might happen that 
the period ends like inflation did, and perhaps the vacuum period is followed by 
a curvature dominated period. If so, the Weyl entropy might again rise to new 
heights, increasing towards a value where all the matter is collected in black 
holes. 
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